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^— ^ ■ These notes provide an introduction to the theory of random matrices. The central quantity 

^y^ ' studied is T{a) = det (1 — K) where K is the integral operator with kernel 

1 sin7r(x — y) 



y 



■xi{y)- 



C_) • Here / = [J, {a2j~i,a2j) and xi{v) is the characteristic function of the set /. In the Gaussian 



Unitary Ensemble (GUE) the probability that no eigenvalues lie in I is equal to T{a). Also T{a) 
is a tau-function and we present a new simplified derivation of the system of nonlinear completely 
integrable equations (the aj 's are the independent variables) that were first derived by Jimbo, Miwa, 
Mori, and Sato in 1980. In the case of a single interval these equations are reducible to a Painleve 
^ . V equation. For large s we give an asymptotic formula for E2{n; s), which is the probability in the 

CO ' GUE that exactly n eigenvalues lie in an interval of length s. 
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I. INTRODUCTION 

These notes provide an introduction to that aspect of the theory of random matrices dealing with the distribution 
of eigenvalues. To first orient the reader, we present in Sec. O some numerical experiments that illustrate some of 



flt"^' the basic aspects of the subject. In Sec. Ill we introduce the invariant measures for the three "circular ensembles" 

(~| , involving unitary matrices. We also define the level spacing distributions and express these distributions in terms of a 

;iV ' particular Fredholm determinant. In Sec. IV we explain how these measures are modified for the orthogonal polynomial 

. !^ " ensembles. In Sec. M we discuss the universality of these level spacing distribution functions in a particular scaling 

S^ limit. The discussion up to this point (with the possible exception of Sec. |^ follows the well-known path pioneered 

5—1 _ by Hua, Wigner, Dyson, Mehta and others who first developed this theory (see, e.g., the reprint volume of Porter |3J] 



and Hua |17|). This, and much more, is discussed in Mehta's book [p5[ — the classic reference in the subject. 

An important development in random matrices was the discovery by Jimbo, Miwa, Mori, and Sato [El| (hereafter 
referred to as JMMS) that the basic Fredholm determinant mentioned above is a r-function in the sense of the Kyoto 
School. Though it has been some twelve years since p^ was published, these results are not widely appreciated by 
the practitioners of random matrices. This is due no doubt to the complexity of their paper. The methods of JMMS 



are methods of discovery; but now that we know the result, simpler proofs can be constructed. In Sec. VI we give 
such a proof of the JMMS equations. Our proof is a simplification and generalization of Mehta's [g^ simplified proof 
of the single interval case. Also our methods build on the earlier work of Its, Izergin, Korepin, and Slavnov ||l^ and 
Dyson ]13] . We include in this section a discussion of the connection between the JMMS equations and the integrable 
Hamiltonian systems that appear in the geometry of quadrics and spectral theory as developed by Moser |^ . This 
section concludes with a discussion of the case of a single interval (viz., probability that exactly n eigenvalues lie in 
a given interval). In this case the JMMS equations can be reduced to a single ordinary differential equation — the 
Painleve V equati on. 

Finally, in Sec. |VH| we discuss the asymptotics in the case of a large single interval of the various level spacing 
distribution functions |H,pq,E8[ . In this analysis both the Painleve representation and new results in Toeplitz/ Wiener- 
Hopf theory are needed to produce these asymptotics. We also give an approach based on the asymptotics of the 
eigenvalues of the basic linear integral operator |l4,G3,Ra] . These results are then compared with the continuum model 
calculations of Dyson Q . 



II. NUMERICAL EXPERIMENTS 

The Gaussian orthogonal ensemble (GOE) consists oi N x N real symmetric matrices whose elements (subject to 
the symmetric constraint) are independent and identically distributed Gaussian random variables of mean zero and 
variance one. One can use a Gaussian random number generator to produce a "typical" such matrix. Given this 
matrix we can diagonalize it to produce our "random eigenvalues." Using the software MATHEMATICA, 25 such 
100 X 100 GOE matrices were generated and Fig. 1 is a histogram of the density of eigenvalues where the x-axis has 
been normalized so that all eigenvalues lie in [—1, 1]. Also shown is the Wigner semicircle law 

Pw{x) = -Vl - x^. (1) 

Given any such distribution (or density) function, one can ask to what extent is it "universal." In Fig. 2 we plot the 
same density histogram except we change the distribution of matrix elements to the uniform distribution on [—1, 1]. 
One sees that the same semicircle law is a good approximation to the density of eigenvalues. See p5[ for further 
discussion of the Wigner semicircle law. 

A fundamental quantity of the theory is the (conditional) probability that given an eigenvalue at a, the next 
eigenvalue lies between b and b + db: p(0; a, b) db. In measuring this quantity it is usually assumed that the system is 
well approximated by a translationally invariant system of constant eigenvalue density. This density is conveniently 
normalized to one. In the translationally invariant system p(0; a, b) will depend only upon the difference s := b — a. 
When there is no chance for confusion, we denote this probability density simply by p{s). Now as Figs. 1 and 2 clearly 
show, the eigenvalue density in the above examples is not constant. However, since we are mainly interested in the 
case of large matrices (and ultimately N — > oo), we take for our data the eigenvalues lying in an interval in which 
the density does not change significantly. For this data we compute a histogram of spacings of eigenvalues. That is 
to say, we order the eigenvalues Ei and compute the level spacings Si := -E^+i — Ei. The scale of both the cc-axis 
and y-axis are fixed once we require that the integrated density is one and the mean spacing is one. Fig. 3 shows the 
resulting histogram for 20, 100 x 100 GOE matrices where the eigenvalues were taken from the middle half of the 
entire eigenvalue spectrum. The important aspect of this data is it shows level repulsion of eigenvalues as indicated 
by the vanishing of p{s) for small s. Also plotted is the Wigner surmise 



Pwis) = -sexp(^--s j (2) 

which for these purposes numerically well approximates the exact result (to be discussed below) in the range < s < 3. 
In Fig. 4 we show the same histogram except now the data are from 50, 100 x 100 real symmetric matrices whose 
elements are iid random variables with uniform distribution on [—1, 1]. 

In computing these histograms, the spacings were computed for each realization of a random matrix, and then 
the level spacings of several experiments were lumped together to form the data. If one first forms the data by 
mixing together the eigenvalues from the random matrices, and then computes the level spacing histogram the results 
are completely different. This is illustrated in Fig. 5 where the resulting histogram is well approximated by the 
Poisson density exp(— s) (see Appendix 2 in (25| for further discussion of the superposition of levels). There are other 
numerical experiments that can be done, and the reader is invited to discover various aspects of random matrix theory 
by devising ones own experiments. 

III. INVARIANT MEASURES AND LEVEL SPACING DISTRIBUTIONS 

A. Preliminary Remarks 

In classical statistical mechanics, the microcanonical ensemble is defined by the measure that assigns equal a priori 
probability to all states of the given system (which in turn is defined by specifying a Hamiltonian on phase space) 
of fixed energy E and volume V. The motivation for this measure is that after specifying the energy of the system, 
every point in phase space lying on the energy surface should be equally likely since we have "no further macroscopic 
information." In the applications of random matrix theory to quantum systems, the Hamiltonian is modeled by a 
matrix H. However, in this case we give up knowledge of the system, i.e. H itself is unknown. Depending upon the 
symmetries of the system, the sets of possible 7?'s are taken to be real symmetric matrices, Hermitian matrices, or 
self-dual Hermitian matrices ("quaternion real") |^,^ 25 1. The question then is what measure do we choose for each 
of these sets of i/'s? 



What we would intuitively like to do is to make each H equally likely to reflect our "total ignorance" of the system. 
Because these spaces of i/'s are noncompact, it is of course impossible to have a probability measure that assigns equal 
a priori probability. It is useful to recall the situation on the real line R. If we confine ourselves to a finite interval 
[a, 6], then the unique translationally invariant probability measure is the normalized Lebesgue measure. Another 
characterization of this probability measure is that it is the unique density that maximizes the information entropy 

S[p] = - p{x) log p{x)dx. (3) 



On R the maximum entropy density subject to the constraints E{1) = 1 and E{x'^) — a^ is the Gaussian density of 
variance ct^ . The Gaussian ensembles of random matrix theory can also be characterized as those measures that have 
maximum information entropy subject to the constraint of a fixed value of E{H*H) p|,|3q|. The well-known explicit 
formulas are given below in Sec. |V[ 

Another approach, first taken by Dyson H] and the one we follow here, is to consider unitary matrices rather than 
hermitian matrices. The advantage here is that the space of unitary matrices is compact and the eigenvalue density 
is constant (translationally invariant distributions). 

B. Haar Measure for U{N) 

We denote by G = U{N) the set of A^ x A^ unitary matrices and recall that dimp^[/(A^) — A^. One can think of 
U{N) as an A^^-dimensional submanifold of R^ under the identification of a complex N x N matrix with a point 
in R^ . The group G acts on itself by cither left or right translations, i.e. iix go £ G then 

Lgo- 9 ^ 909 and Rg^ : g ^ ggo . 

The normalized Haar measure fi2 is the unique probability measure on G that is both left- and right-invariant: 

ti2igE) = ^i2{Eg) ^ fi2{E) (4) 

for all g G G and every measurable set E (the reason for the subscript 2 will become clear below) . Since for compact 
groups a measure that is left-invariant is also right-invariant, we need only construct a left-invariant measure to 
obtain the Haar measure. The invariance (m reflects the precise meaning of "total ignorance" and is the analogue of 
translational invariance picking out the Lebesgue measure. 

To construct the Haar measure we construct the matrix of left-invariant 1-forms 

ng=g-'dg, geG, (5) 

where fig is anti-Hermitian since g is unitary. Choosing A^^ linearly independent 1-forms ujij from fig, we form the 
associated volume form obtained by taking the wedge product of these oj^'s. This volume form on G is left-invariant, 
and hence up to normalization, it is the desired probability measure. 

Another way to contract the Haar measure is to introduce the standard Ricmannian metric on the space oi N x N 
complex matrices Z = (zij): 

N 

(ds)^ = tr{dZdZ*) = J2 M-^yf • 

We now restrict this metric to the submanifold of unitary matrices. A simple computation shows the restricted metric 
is 

{ds)' = tT{ngn;). (6) 

Since fig is left-invariant, so is the metric (ds)"^. If we use the standard Riemannian formulas to construct the volume 
element, we will arrive at the invariant volume form. 

We are interested in the induced probability density on the eigenvalues. This calculation is classical and can be 
found in |17]j37|l . The derivation is particularly clear if we make use of the Riemannian metric (|6|) . To this end we 
write 



g = XDX-\ geG, (7) 

where X is a unitary matrix and I? is a diagonal matrix whose diagonal elements we write as exp(i(pfe). Up to an 
ordering of the angles ipk the matrix D is unique. We can assume that the eigenvalues are distinct since the degenerate 
case has measure zero. Thus X is determined up to a diagonal unitary matrix. If we denote by T{N) the subgroup of 
diagonal unitary matrices, then to each g E G there corresponds a unique pair {X,D), X E G/T{N) and D E T{N). 
The Haar measure induces a measure on G/T{N) (via the natural projection map). Since 

X*dgX = VLxD - DVtx + dD , 

we have 

(ds)2 = tr {ng^l) = tr (dgdg*) = tr {X*dgXX*dg* X) 
= tr ([J7x, £>] [fix, -D]*) + tr {dD dD*) 

N N 

= ^ |(5Xfcf (exp(z(/7fc) ~ Gx^[iipt))f + y^(dy^fc)^ 

k,l=\ fc=l 

where fix — (6Xki). Note that the diagonal elements iX^^ do not appear in the metric. Using the Ricmannian 
volume formula we obtain the volume form 

^^3=^-^11 I^^P(*'f^i) ~ exp(i(pfc)|^ d^x ■ ■ -dipN (8) 

j<k 

where ujx = const J|-^j, SXjk- 

We now integrate over the entire group G subject to the condition that the elements have their angles between ip^ 
and ipk + dipk to obtain 

Theorem 1 The volume of that part of the unitary group U{N) whose elements have their angles between ip^. and 
ipk + dipk is given by 

-PAf2(t/5i,. •■,</5Ar)d<^i ■ • ■dfySA' = Cjv2 JJ |exp(i(pj) -exp(i(pfc)| dipi---d(pN (9) 

j<k 

where Gn2 is a normalization constant. 

C. Orthogonal and Symplectic Ensembles 

Dyson |10| , in a careful analysis of the implications of time-reversal invariance for physical systems, showed that 
(a) systems having time-reversal invariance and rotational symmetry or having time-reversal invariance and integral 
spin are characterized by symmetric unitary matrices; and (b) those systems having time-reversal invariance and 
half-integral spin but no rotational symmetry are characterized by self-dual unitary matrices (see also Chp. 9 in p5[). 
For systems without time-reversal invariance there is no restriction on the unitary matrices. These three sets of 
unitary matrices along with their respective invariant measures, which we denote by Efs{N), (3 = 1,4,2, respectively, 
constitute the circular ensembles. We denote the invariant measures by /i^, e.g. /X2 is the normalized Haar measure 



discussed in Sec. Ill B 



A symmetric unitary matrix S can be written as 

S = V^V, V E U{N). 

Such a decomposition is not unique since 

V ^ RV, ReO{N), 

leaves S unchanged. Thus the space of symmetric unitary matrices can be identified with the coset space 

U{N)/0{N) . 

The group G = U{N) acts on the coset space U{N)/0{N), and we want the invariant measure ^i. If tt denotes the 
natural projection map G -^ G/H, then the measure fii is just the induced measure: 



MB) ^ f,2{7T-\B)) , 

and hence Ei{N) can be identified with the pair {U{N)/0{N), ni). 

The space of self-dual unitary matrices can be identified (for even N) with the coset space 

U{N)/Sp{N/2) 

where Sp{N) is the symplectic group. Similarly, the circular ensemble E4^{N) can be identified with 
{U{N)/ Sp{ N/2), f i4) where 114 is the induced measure. 



As in Sec. [II B we want the the probability density Pni3 on the eigenvalue angles that results from the measure ug 
ioT P — 1 and P — 4. This calculation is somewhat more involved and we refer the reader to the original literature PJ17|] 
or to Chp. 9 in |23]. The basic result is 

Theorem 2 In the ensemble Ep{N) ([3 — 1, 2, A) the probability of finding the eigenvalues exp(i(y9j) of S with an angle 
in each of the intervals {0j, 6j + dOj) (j — I, . . . , N) is given by 

PNpiOu...,0N)dei---deN^CNp n \eMiOi)-exp{i0j)f dei---deN (lo) 

i<e<j<N 
where C^fj is a normalization constant. 
The normalization constant follows from 
Theorem 3 Define for N eN and /3 G C 

^n{P) - (27r)-^ / • • • / TT |exp(z0,) - exp{tek)f d9i ■■■dON (11) 

then 

{T{l + fll2)f 

The integral ([ll| ) has an interesting history, and it is now understood to be a special case of Selberg's integral (see, 
e.g., the discussion in p5|]). 

D. Physical Interpretation of the Probabihty Density Pjv/3 
The 2D Coulomb potential for N unit like charges on a circle of radius one is 

WN{9i,...,eN)^- Yl log|exp(z0,)-exp(z0fe)| . (13) 

i<i<fc<JV 

Thus the (positional) equilibrium Gibbs measure at inverse temperature < /? < 00 for this system of N charges with 
energy Wn is 

exp(-/3W^^(gi,...,0Ar)) 

^n{P) ■ ^ ^ 

For the special cases of /? — 1,2,4 this Gibbs measure is the probability density of Theorem 0. Thus in this math- 
ematically equivalent description, the term "level repulsion" takes on the physical meaning of repulsion of charges. 
This Coulomb gas description is due to Dyson 0, and it suggests various physically natural approximations that 
would not otherwise be so clear. 



E. Level Spacing Distribution Functions 

For large matrices there is too much information in the probabihty densities Piqp{di, . . . ^Oiq) to be useful in 
comparison with data. Thus we want to integrate out some of this information. Since the probability density 
PNp{di, ■ ■ • 7 On) is a completely symmetric function of its arguments, it is natural to introduce the n-point correlation 
functions 

N\ f"^^ Z"^'" 

Rnpi^l,- ■ ■ T^n) — JJZ TT / •■■/ PnI3{0i,. . . ,6iS[)d9n^i- ■ -dOisi . (15) 

{1\ — n). Jq Jq 

The case f3 — 2 is significantly simpler to handle and we limit our discussion here to this case. 



1 / ^ 

-det{KN(9,,0k) 

Lemma 1 where 



PiV2(0i,...,0jv) = ^det(if^(0j,0fe) _ ) (16) 



1 sin(iV(g,--gfc)/2) 
^"('-''=)=2;^sin((^,-^,)/2) ■ (^') 

Proof: 

Recalling the Vandermonde determinant we can write 

Jl |exp(i6'j) - exp(i6'fe)|^ = det(M'^) det(M) (18) 

]<k 

where Mjk = exp (i{j ~ l)Ok). A simple calculation shows 

{M^M)^^ = 27:D,,KNie,,0k)Dkk (19) 

where D is the diagonal matrix with entries exp(z(Af — 1)9 j/2). Except for the normalization constant the lemma 
now follows. Getting the correct normalization constant requires a little more work (see, e.g., Chp. 5 in 123]). ■ 

From this lemma and the combinatoric Theorem 5.2.1 of [E5| follows 



Theorem 4 Let i?n2(^i, ■ • ■ , 0n) be the n-point Junction defined by Al3) for the circular ensemble E2{N); then 

Rn2{0i, . . . ,0„) = det (KN{9j,6k)'\\ (20) 

where KN{9j,9k) is given by Jllij). 

We now discuss the behavior for large N. The 1-point correlation function 

RiAOi) = ^ (21) 

is just the density, p, of eigenvalues with mean spacing D = 1/ p. As the size of the matrices goes to infinity so does 
the density. Thus if we are to construct a meaningful limit as A^ ^ oo wc must scale the angles 9j . This motivates 
the definition of the scaling limit 

p -^ oo, 9j -^ 0, such that Xj := p9j G R is fixed. (22) 

We will abbreviate this scaling limit by simply writing N —^ oo. In this limit 

Rn2{xi,...,Xn)dxi---dXn-^ lim i?„ (6*1 ,..., 6I„) ^6*1 ••• (i6'„ (23) 

where we used the slightly confusing notation of denoting the scaling limit of the n-point functions by the same 
symbol. From Theorem follows 



Theorem 5 In the scaling limit l[2^) the n-point functions become 

R.n2{xi,...,Xn) = A.etiK{Xj,X},) ] (24) 

where the kernel K{x,y) is given by 



j,k=i 



1 sinTrfx — v) , , 

TT X — y 

The three sets of correlation functions 

f/3:={i?«/3(a;i,...,a:„);x, eRj^^i, /3= 1,2,4, (26) 

define three different statistics (called the orthogonal ensemble, unitary ensemble, symplectic ensemble, respectively) 
of an infinite sequence of eigenvalues (or as we sometimes say, levels) on the real line. 

We now have the necessary machinary to discuss the level spacing correlation functions in the ensemble £2- We 
denote by X the union of m disjoint sub-intervals of the unit circle: 

X = XiU---UX™. (27) 

We begin with the probability of finding exactly ni eigenvalues in interval Xi,. . . , nm eigenvalues in interval 2m in 
the ensemble E2{N). We denote this probability by E2N{ni, ■ ■ ■ , Um',!) and we will let A^ — > 00 at the end. 
If XA denotes the characteristic function of the set A and n :~ ui + ■ ■ ■ + Um, then the probability we want is 



£'2Ar(ni, . . . ,nm;I) = ,, / dOi ■ ■ ■ d0N PN2{di, ■ ■ ■ ,0n) 

\ni ■■■nmN -nj Jq Jq 



rii ni+n2 niH hrim 

X n ^^1 ^^n ) n ^^^ (^n ) • • • n ^^^ (^^'" ) 

N 

X n (1 - ^^(^^■)) • (28) 

j=n+l 



We define the quantities 



<.27r 
■ / dOn Rn{Oi, . . . ,9n) 
JQ 



2tt 

d9i 



"1 ni+n2 »ilH h"™ 

X n ^^1 (^n ) n ^^^ (^^-^ ) • • • n ^^™ ^^^^ ) ■ (^9) 

Jl=l J2=ni+1 j„=niH hi 



The idea is to expand the last product term in (E8|) involving the characteristic function xx and to regroup the terms 
according to the number of times a factor of xi appears. Doing this one can then integrate out those 9k variables that 
do not appear as arguments of any of the characteristic functions and express the result in terms of the quantities 
rni---nm- To recognize the resulting terms we define the Fredholm determinant 

1-J2^jKn{9,9')xx,{0') 

where "J2T=i ^j-^n{9,9')xIj{9')" means the operator with that kernel and A is the ?7i-tuple (Ai, . . . , Am)- A slight 
rewriting of the Fredholm expansion gives 

Jl H \-JTn=J 

The expansion above is then recognized to be proportional to 

d'WN{T,X) 



uA^ • • • OArn 

The scaling limit N ^ oo can be taken with the result: 



Ai=-=A,„ = 1 



Theorem 6 Given m disjoint open intervals Ik — (a2fe-i,a2fe) C R, let 

/:=/iU---U/,„. 



(30) 



The probability E2{ni, . . . , rim] I) in the ensemble £2 that exactly n^ levels occur in interval Ik (7c = 1, ... , m) is given 
by 



£'2(^1,.. .,nm,I) 



(-1)" d"D{I;X) 



Til! •••«„! aA"i •••9A'i 



Ai=-=A,„ = 1 



(31) 



whe 



D{I; A) - dot I 1 - ^ X,K{x, y)xi, (y) | , 



(32) 



K{x, y) is given by (pq^, and rt := ni + • • • + n„ 



In the case of a single interval / = (a, 6), we write the probability in ensemble £p of exactly n eigenvalues in / as 
Ep{n; s) where s := b ~ a. Mehta |23,Effl has shown that if we define 



D±{s;X) =det{l-XK±) 
where K± are the operators with kernels K{x, y) ± K{—x, y), and 

(-1)" d"D±{s;X) 



then£;i(0;s) =£+(0;s), 



E±{n;s) = 



i\ a A" 



A=l 



E+{n;s) ^ Ei{2n; s) + Eii2n - l;s), n > 



(33) 



(34) 



(35) 



and 



E_{n;s) = Ei{2n; s) + Eii2n + 1; s), n> , 



Ei{n\s) = -{E+{n;2s) + E^{n;2s)) , n>0 . 



(36) 



(37) 



Using the Fredholm expansion, small s expansions can be found for ii^^(n; s). We quote ||25| , P6| here only the results 
for n = 0: 



Ei{0;s) = 1-s 



TT^S^ TT^'s^ 



36 1200 



+ 0(56) 



36 675 



- + 0(.«) 



^^(«^^) = l-^ + ^02^ + «(^)- 



(38) 



The conditional probability in the ensemble £(3 of an eigenvalue between b and b + db given an eigenvalue at a, is 
given [Q by P/3(0; s) ds where 



P/3(0;s) 



d^Ep{0;s) 
ds^ 



(39) 



Using this formula and the expansions ( pq ) we see that pg{Q;s) — 0{s^), making connection with the numerical 
results discussed in Sec. ||. Note that the Wigner surmise (g) gives for small s the correct power of s for pi{0; s), but 
the slope is incorrect. 



IV. ORTHOGONAL POLYNOMIAL ENSEMBLES 

Orthogonal polynomial ensembles have been studied since the 1960's, e^. [OLbut recently interest has revived 
because of their application to the matrix models of 2D quantum gravity n,@,n3.n6[ . Here we give the main results 
that generalize the previous sections. 

The orthogonal polynomial ensemble associated to V assigns a probability measure on the space oi N x N hermitian 
matrices proportional to 

exp (-Tr {V{M))) dM (40) 

where V{x) is a real- valued function such that 

wv{x) — exp(— F(x)) 

defines a weight function in the sense of orthogonal polynomial theory. The quantity dM denotes the product of 
Lebesgue measures over the independent elements of the hermitian matrix M . Since Tr {V{M)) depends only upon 
the eigenvalues of M, we may diagonalize M 

M = XDX* , 

and as before integrate over the "X" part of the measure to obtain a probability measure on the eigenvalues. Doing 
this gives the density 

/ N \ 

Pn{xi,...,xn)= W {xk-xtfex:p\-^V{xj)\. (41) 

l<k<e<N \ j=l J 

If we introduce the orthogonal polynomials 

Pmix)Pn{x)wvix)dx = 6mn, m, 71 = 0,1,... (42) 



/R 

and associated functions 



(Pm{x) = e:>q){-V{x)/2)pmix) , 
then the probability density ( [4l| ) becomes 

Pn{xi,...,xn) = ^ (^det((^j_i(xfc))|^^^J 



— detiKN{x„Xk)) (43) 



where 



JV-l 

Kn{x,v) = ^ Lpj{x)Lpj{y) 
j=o 

_ kN-l ipNix)ipN-l{y) - (pN-lix)(pNiy) 

kpf x-y 



for X ^ y 



= — ((p'pf{x)ipN{x)-ip'pf^i{x)ipjy{x)) for x = y. (44) 

The last two equalities follow from the Christoffel-Darboux formula and the fc„ are defined by 

Pn\X) = KjiX + • • • , fc„ > U. 

Using the orthonormality of the tpj {x) 's one shows exactly as in |2^ that 

R„{xi,...,Xn) ■■ ^ jrrz ^-7 ■■ Pn{xi, ■ ■ ■ ,Xn) dx^+l ■ ■ ■ dxjs 

= det(^KN{xj,Xk)\2k^^J ■ (45) 



In particular, the density of eigenvalues is given by 

Pat (a;) = Kn{x,x) . (46) 

Arguing as before, we find the probability that an interval / contains no eigenvalues is given by the determinant 

det{l-KN) (47) 

where Kn denotes the operator with kernel 

KN{x,y)xi{y) 

and Kpf{x, y) given by ([44|). Analogous formulas hold for the probability of exactly n eigenvalues in an interval. We 
remark that the size of the matrix N has been kept fixed throughout this discusion. 

V. UNIVERSALITY 

We now consider the limit as the size of the matrices tends to infinity in the orthogonal polynomial ensembles of 



Sec. IV. Recall that we defined the scaling limit by introducing new variables such that in these scaled variables 



the mean spacing of eigenvalues was unity. In Sec. Ill the system for finite A^ was translationally invariant and had 
constant mean density N/2n. The orthogonal polynomial ensembles are not translationally invariant (recall (^6|) so 
we now take a fixed point xq in the support of pjv (x) and examine the local statistics of the eigenvalues in some small 
neighborhood of this point xq. Precisely, the scaling limit in the orthogonal polynomial ensemble that we consider is 

A^ — > oo, Xj — > Xq such that ^j := pj^{xo){xj — a;o) is fixed. (48) 

The problem is to compute Kn{x, y) dy in this scaling limit. From ( [44| ) one sees this is a question of asymptotics of 
the associated orthogonal polynomials. For weight functions wyix) corresponding to classical orthogonal polynomials 
such asymptotic formulas are well known and using these it has been shown |13| , p2| that 

J^ f ^w 1 sin7r(^-^0 , 

Kn[x, y)dy^ — d£, . (49) 

TT £,-£,' 

Note this result is independent of xq and any of the parameters that might appear in the weight function wv{x). 
Moore [Q has given heuristic semiclassical arguments that show that we can expect (|9|) in great generality. 

There is a growing literature on asymptotics of orthogonal polynomials when the weight function wy {x) has poly- 
nomial V , see ^^ and references therein. For example, for the case of V{x) — x^ Nevai [Q has given rigorous 
asymptotic formulas for the associated polynomials. It is straightforward to use Nevai's formulas to verify that (|49| ) 
holds in this non-classical case. 

There are places where (|49| ) will fail and these will correspond to choosing the a;o at the "edge of the spectrum" 
P,pO[ (in these cases xq varies with NY This will correspond to the double scaling limit discovered in the matrix 
models of 2D quantum gravity |7|,|8|, |l5| , [16| . Different multicritical points will have different limiting Kn{x, y) dy [|6|,pO| . 



VI. JIMBO-MIWA-MORI-SATO EQUATIONS 

A. Definitions and Lemmas 

In this section we denote by a the 2TO-tuple (ai, . . . , a2m) where the aj are the endpoints of the intervals given in 
Theorem H and by da exterior differentiation with respect to the Uj {j — 1, . . . , 2rn). We make the specialization 

Aj = A for j 

which is the case consider in [gl|. It is not difficult to extend the considerations of this section to the general case. 
We denote by K the operator that has kernel 

XK{x,y)xi{y) (50) 
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where K{x, y) is given by (pq). It is convenient to write 



\K{x, y) = ^— (51) 



where 



A(x) = sinTra;. 

The operator K acts on L^(R), but can be restricted to act on a dense subset of smooth functions. From calculus we 
get the formula 

K=i-iyKix,a,)Siy-a,) (52) 

oaj 

where S{x) is the Dirac delta function. Note that in the right hand side of the above equation we are using the shorthand 
notation that "A{x, y)" means the operator with that kernel. We will continue to use this notation throughout this 
section. 

We introduce the functions 

Q{x; a) = (1 - K)-^Aix) = f p{x, y)A{y) dy (53) 

and 

P{x- a) - (1 - Kr^A\x) - / p{x, y)^{y) dy (54) 

JR 

where p{x, y) denotes the distributional kernel of (1 — K)^^ . We will sometimes abbreviate these to Q{x) and P{x), 
respectively. It is also convenient to introduce the resolvent kernel 

i? = (1 ~ K)-^K . 

In terms of kernels, these are related by 

p{x,y) = S{x~y) + R{x,y). 

We define the fundamental 1-form 

uj{a) := da log D [I ■,X). (55) 

Since the integral operator K is trace-class and depends smoothly on the parameters a, we have the well known result 

Lu{a) = -Tr ((1 - K)-^daK) . (56) 

Using ( p2| ) this last trace can be expressed in terms of the resolvent kernel: 

2m 

u^{a) = -Y.{-iyR{a„a,)da, (57) 

j=l 

which shows the importance of the quantities R{aj, Uj). A short calculation establishes 

d 



daj 



(1 - K)-^ = {-lYR{x, aj)p{aj,y) . (58) 



We will need two commutators which we state as lemmas. 
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Lemma 2 If D = -^ denotes the differentiation operator, then 



2m 

[D, (1 - Xy'] = - ^ (-!)■'■ i?(a;, a,>(a„ y) . 



Proof: 
Since 



[D, (1 - K)-'] ={1- K)-' [D, K] (1 - K)-' , 
we begin by computing [D, K]. An integration by parts shows that 

[Z?,if] = -^(-l)-'"i^(a;,a,)<5(y-a,) 



where we used the property 

dK{x,y) dK{x,y) 
dx dy 



= 



satisfied by our K{x,y) and the well known formula for the derivative of x/(a;). The lemma now follows from the fact 
that 

Ri^iV) ^ / p{x,z)K{z,y)dz = / K{x, z)p{z,y) dz . 

JR JR 



Lemma 3 // M^ denotes the multiplication by x operator, then 

[M,,{1^K)-'] =Q{x) {l^KY'A'xi{y)-P{x) {l~Ky^Axi{y) 
where K* denotes the transpose of K , and also 

M,,il-K)-'] ={x-y)R{x,y). 



Proof: 
We have 



[M,, K] - {A{x)A'{y) - A'{x)A{y)) xM (1 - KY 



From this last equation the first part of the lemma follows using the definitions of Q and P. The alternative expression 
for the commutator follows directly from the definition of p(x, y) and its relationship to R(x, y). ■ 

This lemma leads to a representation for the kernel R{x,y) |lq]: 



Lemma 4 If R{x,y) is the resolvent kernel of ([5jj and Q and P are defined by ( |5q j and m^, respectively, then for 
x,y Cz I we have 

_ Q(x; a)P{y; a) - P{x; a)Q{y; a) 

j^[Xj y) — 1 X r y 1 

x-y 



and 



R{x, x) = —;—{x; a) P{x; a) ;— (a;; a) Q{x; a) . 

dx dx 
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Proof: 

Since K{x, y) — K{y, x) we have, on /, 

(1 - K'r'Axi = (1 - K)-'Axi = (1 - K)-'A 

(the last since the kernel of K vanishes for y ^ I). Thus (1 — K*^)^^Axi = Q on I, and similarly (1 — K*)~^A'xi — P 
on /. The first part of the lemma then follows from Lemma 0. The expression for the diagonal follows from Taylor's 
theorem. ■ 



We remark that Lemma W used only the property that the kernel K{x, y) can be written as ( |5l| ) and not the specific 
form of A{x). Such kernels arc called "completely integrable integral operators" by Its et. al. ||l^ and Lemma || is 
central to their work. 

B. Derivation of the JMMS Equations 

We set 

qj = qj{a) = lim Q{x]a) and pj = pj{a) — lim P{x;a), j = l,...,2m. (59) 



Specializing Lemma we obtain immediately 



Referring to (|8|) we easily deduce that 

dq. 



R{a.j,ak)^- — , J^k- (60) 



-^ = {-lfR{a,,ak)qk. J^k, (61) 

oak 



and 






i~l)''R{a,,ak)pk, J^k. (62) 



Now 



^ ^ D (1 - K)-^ A{x) 
ax 



Thus 



Similarly, 



= (1 - Ky^ DA{x) + \d, (1 - Ky^] A{x) 
= (1 - A^-^ A'{x) - Y,{-iyR{x, ak)qk ■ 



— {aj;a) ^ p^ ~^{-iyR{aj,ak)qk. (63) 



riP r n 

= (l-Ky^A"{x)+\D,{l-Ky^ A'{x) 
ax I i 

= -7^2 (1 _ Ky^ A{x) - ^(-l)'=i?(x, ak)pk 
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Thus 



dP_ 

dx 



{aj;a) = -Tr^g^ - ^{-l)''R{aj,ak)pk 



Using ( p3| ) and ( |64| ) in the expression for the diagonal of R in Lemma ^, we find 

R{aj,a.j) = 7r^g| +p^ + ^(-l)''ii(aj, afe)i?(afe, aj)(aj - Ofc) . 



Using 



(|6l|) and (B3) we obtain 



Similarly, 



9(7j dQ{x;a) 



dx 



dQ{x;a) 



do 



daj 



Pj - ^{-'^)''R{aj,ak)qk 



k^J 



dp^ 






da 



(64) 



(65) 



(66) 



(67) 



Equations (60)-(p3) and (|65|)-(|67|) are the JMMS equations. We remark that they appear in slightly different form 
in |2l|] due to the use of sines and cosines rather than exponentials in the definitions of Q and P. 



C. Hamiltonian Structure of the JMMS Equations 

To facilitate comparison with |^,0 we introduce 

_ i _l 

<l2j - ~-X2j , q2j+l - -X2J + 1 , 

P2j = -iy2j , P2J + 1 = y2j + l , 

Mjk := -^(xjyk -XkVj), 

fc=i 

In this notation, 



m: 



jk 



(68) 



^(a) = ^Gj{x,y)daj. 
j 

If we introduce the canonical symplectic structure 

{xj,xk} = {yo,yk} = 0, {xj^yk} = Sjk , (69) 

then as shown in Moser pit] 

Theorem 7 The integrals Gj{x,y) are in involution; that is, if we define the symplectic structure by MS) we have 

{Gj , Gfe} = for all j,k = 1,..., 2m. 

Furthermore as can be easily verified, the JMMS equations take the following form pll]: 
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Theorem 8 // we define the Hamiltonian 



2m 

Lj{a) — y. Gj{x,y)daj , 



then Eqs. mi), n62), fl6q) and (61) are equivalent to Hamilton's equations 

daXj = {xj,u{a)} and daijj = {yj,uj{a)} . 

In words, the flow of tlie point (a;, y) in the "time variable" aj is given by Hamilton's equations with Hamiltonian Gj. 
The (Frobenius) complete integrability of the JMMS equations follows immediately from Theorems and ||. We 
must show 

da{xj,uj} — and da{yj,u>} — 0. 

Now 

da {Xj , U} = {daXj , Uj} + {Xj , daUj} , 

but daW = since the Gfc's are in involution. And we have 

{daXj,uj} = ^ {{{xj,Gk},Gi} - {{xj,Gi},Gk})dak A dag 

k<f. 

which is seen to be zero from Jacobi's identity and the involutive property of the Gj's. 



D. Reduction to Painleve V in the One Interval Case 

1. The (j{x;\) differential equation 

We consider the case of one interval: 

m — 1, ai = —t, 02 = t, with s := 2t . (70) 

Since p{x, y) is both symmetric and even for x,y £ I, we have 92 = ^<li and p2 = Pi- Introducing the quantity 

^1 — / p{^t,x)exp{—iTrx)dx, 



we write 



qi = — -(ri -ri) and Pi = -r- (^i + ri ) . 
Zm 2 



Specializing the results of Sec. VI B to m = 1 we have 



Lo{a)^-2R{t,t)dt, (71) 

R{-t,t)^-p,p,^^^{rl-rj), (72) 

dqi dqi dqi , ouf ^ +^ 

at oai oa2 

^=^2(7i+2i?(-i,t)pi, 

^=z7rri+2i?(~i,t)ri, (73) 

R(t, t) = TT^q? +pI- 2tR{-t, tf 

= Xnn + ^^{rl-riy. (74) 

15 



A straightforward computatfon from (72)-(|74|) shows 

d 



f^{tRi-t,t)) ^X^irl), (75) 

l{tR{t,t))=X\n\\ (76) 

j^Rit,t)^2iR{-t,t)f . (77) 



Eq. (frn) is known as Gaudin's relation and Eqs. (|75[) and (iTq) are identities derived by Mehta p7| (see also Dyson [h2| 



^^^l^d.) , (79) 



in his proof of the one interval JMMS equations. Here we made the JMMS equations central and derived (75)~(|77|) 
as consequences. 

These equations make it easy to derive a differential equation for 

a(x;X) := -2tR{t,t) = x--\ogD(-:X), where x = 2TTt. (78) 

dx n' 

We start with the identity 

kl^=(5R(r?))%(3(r?))^ 

and define temporarily a{t) := tR{t,t) and h{t) := tR{~t,t); then ([7^), (fTq), and (|7q ) imply 

Using ( [77|) and its derivative to eliminate b and db/dt, we get an equation for a(i) and hence (t{x; A): 
Theorem 9 In the case of a single interval I — {—t,t) with s = 2t, the Fredholm determinant 

D{s;X) ^det{l-XK) 
is given by 

D{s; X) = exp (J 

where (j{x; X) satisfies the differential equation 

{xa"f + 4 {xa' - a) {xa' -a+ [a'f) = (80) 

with boundary condition as x —^ 

(T(a;;A) = --x-(-)V . (81) 

TT TT 

Proof: Only ( |8l| ) needs explanation. The small x expansion of a{x; A) is fixed from the small s expansion of D{s; A) 
which can be computed from the Neumann expansion. ■ 

The differential equation (^^ is the "cr representation" of the Painleve V equation. This is discussed in |21|] , and in 
more detail in Appendix C of [Q. In terms of the monodromy parameters 9i {i — 0,1, oo) of J2^], ( |80| ) corresponds 
to 00 = 9i = 6ao = which is the case of no local monodromy. For an introduction to Painleve functions see jl^,^ . 

2. The G±{x;X) equations 

Recalling the discussion following Eq. (p3|), we see we need the determinants D±{s; X) to compute Epin; s) for /? = 1 
and 4. Let R± denote the resolvent kernels for the operators 

16 



where {J f){x) — f{—x) and the last equahty of the above equation follows from the evenness of K. Thus 

R± ■.= {l-K±)-'K±^^il±J)R, 



which in terms of kernels is 



Thus 



R±{x,y) ^ ^{R{x,y) ± Ri-x,y)) . 



{R^{t,t) - R+{t,t)f = (R{-t,t)f = ~R{t:t) 



Introducing the analogue of (t{x; A), i.e. 



the above equation becomes 



d X 

(j±{x;X) -.^ x—logD±{-;X), 
ax IT 



cr_ (a;; A) — f7+(x; A) \ d a{x]X) 



(82) 



dx X 
Of course, (t±{x] A) also satisfy 

(j+(x] A) + (T_(a;; A) = (j(x] A) . (83) 



Using ( |82| ) and ( |83| ) and integrating a±{x] X)/x we obtain (the square root sign ambiguity can be fixed from small x 
expansions) 



Theorem 10 Let D±(s;X) be the Fredholm determinants defined by (ujj, then 



1 I f d^ 

logD±{s;\) = -\ogD{s;X)±-J ^ -—\ogD{x; X) dx . (84) 



VII. ASYMPTOTICS 
A. Asymptotics via the Painleve V Representation 

In this section we explain how one derives asymptotic formulas for Efj{n; s) as s — > cxd starting with the Painleve V 
representations of Theorems g and |o[ This section follows Basor, Tracy and Widom Q (see also |2^). We remark 
that the asymptotics of Ep{0;s) as s — > oo was first derived by Dyson ||ll| by a clever use of inverse scattering 
methods. 

Referring to Theorems |^ and ^ one sees that the basic problem from the differential equation point of view is to 
derive large x expansions for 

ao(x):=a(x;l) (85) 

^«(x):=|^(x;l), n=l,2,... (86) 

a±,nix):=^^ix;l), n=l,2,.... (87) 

We point out the sensitivity of these results to the parameter A being set to one. This dependence is best discussed 
in terms of the differential equation ( pO[ ) where it is an instance of the general problem of connection formulas, see 
e.g. [ p2[ and references therein. In this context the problem is: given the small x boundary condition, find asymptotic 
formulas as a; ^ oo where all constants not determined by a local analysis at oo are given as functions of the parameter 
A. If we assume an asymptotic solution for large x of the form a{x) ^ ax^, then (pQ) implies either p = 1 or 2 and 
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a p = 2 then necessarily a = —j. The connection problem for (80) has been studied by McCoy and Tang |gj] who 
show that for < A < 1 one has 

a{x, A) = a{X)x + 6(A) + o(l) 

as a; ^ oo with 

a(A) = - log(l - A) and 6(A) == la^{X) • 

Since these formulas make no sense at A = 1, it is reasonable to guess that 

a{x;l)^-^x\ (88) 

For a rigorous proof of this fact see pS]. It should be noted that in Dyson's work he too "guesses" this leading 
behavior to get his asymptotics to work (this leading behavior is not unexpected from the continuum Coulomb gas 
approximation). Given (p3), and only this, it is a simple matter using (pG) to compute recursively the correction 
terms to this leading asymptotic behavior: 

ao(.)=-i.x^4 + E|^. x^oo (89) 

n— 1 

(c2 — —\,Ci = — |, etc.). Using ( |89| ) in ( [79| ) and ( p^ ) one can efficiently generate the large s expansions for D{s; 1) and 
D±(s; 1) except for overall multiplicative constants. In this instance other methods fix these constants (see discussion 
in HjM). In general this overall multiplicative constant in a r- function is quite difficult to determine (for an example 
of such a determination see |0,^). We record here the result: 

logi?±(s;l)--^7r2s2Tj^s-ilog7rs±ilog2+^log2 + ^C'(-l) + o(l) (90) 

as s ^ cxD where C is the Riemann zeta function. We mention that for < A < 1 the asymptotics of D(s; A) as s — > oo 
are known &p5|. 

One method to determine the asymptotics of cr„(x) as a; — > oo is to examine the variational equations of (|8Q), 
i.e. simply differentiate (|80| ) with respect to A and then set A = 1. These linear differential equations can be solved 
asymptotically given the asymptotic solution (p9) . In carrying this out one finds there are two undetermined constants 
corresponding to the two linearly independent solutions to the first variational equation of (^y). One constant does 
not affect the asymptotics of ai{x) (assuming the other is nonzero!). Determining these constants is part of the 
general connection problem and it has not been solved in the context of differential equations. In [[f[|3q| Toeplitz and 
Wiener-Hopf methods are employed to fix these constants. The Toeplitz arguments of Ej depend upon some unproved 
assumptions about scaling limits, but the considerations of |3q| are completely rigorous and we deduce the following 
result Q for cr„(a;) for all n = 3, 4, . . .: 



n\ exp(nx) 
(^n{x) = - 



l + i(7n-4)i + ^(7n^ + 12n-16)-l + 0(-l) 



(237r)«/2 W2-1 [" ' 8'"" ^'x ' 128'"" ' ^"'" ^^' x^ ' ^'x^'\ ^^^' 

as X —^ cxD. For n = 1,2 the above is correct for the leading behavior but for n = 1 the correction terms have 
coefficients | and r^, respectively, and for n — 2 the above formula gives the coefficient for 1/x but the coefficient 
for 1/x^ is ||. See [Q for asymptotic formulas for <7±^n{x)- 

Since the asymptotics of Efj{0; s) are known, it is convenient to introduce 

Ep{n;s) 

Here we restrict our discussion io j3 — 2 (see Q for other cases). Using (£l|) in ( |3l| ) one discovers that there is a 
great deal of cancellation in the terms which go into the asymptotics of r2(n; s). To prove a result for all n g N by 
this method we must handle all the correction terms in (|9l|) — this was not done in [Q and so the following result was 
proved only for 1 < n < 10: 



exp(n7rs) 
r2(n;s) = S2,„- 



^ n^ 11' 

(2^' + 7) - + TTrA^n' + 48^^ + 229) ^^ + O(-) 
TTs 128 (tts)^ S'^ 



(92) 



iV2 

where 

B2^n = 2-"'-"/2^-("'+")/2 {n - 1)! (n - 2)! • • • 2! l! . 

In the next section we derive the leading term of ( |9^ ) for all n e N. Asymptotic formulas for r^(n; s) {(i — 1,4,±) 
can be found in Q . 
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B. Asymptotics of r2 (n;s) from Asymptotics of Eigenvalues 



The asymptotic formula (92) can also be derived by a completely different method (as was briefly indicated in M). 



If we denote the eigenvalues of the integral operator /C by Ao > Ai > • • • > 0, then 

oo 

det{l-XK) = J|(l- AA,), 
and so it follows immediately from (|3^) that 

(This is formula (5.4.30) in p^.) Thus the asymptotics of the eigenvalues A^ as s ^ oo can be expected to give 
information on the asymptotics of r2(ri; s) as s — > oo. 

It is a remarkable fact that the integral operator K, acting on the interval {—t,t), commutes with the differential 
operator C defined by 

the boundary condition here is that / be continuous at ±t. Thus the integral operator and the differential operator 
have precisely the same eigenfunctions — the so-called prolate spheroidal wave functions, see e.g. p9[. Now Fuchs [14[, 
by an application of the WKB method to the differential equation, and using a connection between the eigenvalues 
Xi and the values of the normalized eigenfunctions at the end-points, derived the asymptotic formula 

1 _ A, ^ ^»+i22'+3/2s»+i/2e— /i! (95) 

valid for fixed i as s ^ oo. Further terms of the asymptotic expansion for the ratio of the two sides were obtained by 
Slepian [||. 

If one looks at the asymptotics of the individual terms on the right side of (^3|) , then we see from ( |95| ) that they all 
have the exponential factor e"'^* and that the powers of s that occur are 

g-n/2-(ji+.-- + i„) ^ 

Thus the term corresponding to ii = 0, i2 = 1, ■■■, in — n — 1 dominates each of the others. In fact we claim this 
term dominates the sum of all the others, and so 



r2 



(n; s) - 1! 2! • • • (n - 1)1 ^-"("+i)/22-»'-"/25-nV2g«^. ^ (96) 



in agreement with (P3). 

To prove this claim, we write 



A,o • • • A,o _ — ./ \ ... A 

^^^"=^^ = (i-A,o)...(i-A,o)+>^(r 



(l-A,o)...(l-A,o) ^ (1-A,J...(1-A,J 

where ii = 0, i^ — 1, . . . , i'^ — n — 1 and the sum here is taken over all (ii, . . . , i„) 7^ (i?, . . . , i*^) with ii < • • • < i„ 
We have to show that 



S'n 



K 



(1-A,J---(1-A,J/ s»V2 

as s — > 00 and we know that this would be true if the sum were replaced by any summand. Write ^ — '^1 +X]2 
where in J^i '^^ have ij < N for all j {N to be determined later) and ^2 i^ the rest of ^ . Since J^i i^ ^ finite sum 
we have 

so we need consider only ^2- I^ ^^Y summand of this we have ij > N for some j, and so for this j 
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1 1 

< 



1 — Ajj 1 — Aat 
and by ( |95| ) this is at most 

„ ^-Af-l/2 „7i-s 

a^v s ' e 

for some constant cat. The product of aU other factors 

1 



1-A, 
appearing in this sumniand is at most 



1 \ n— 1 



,1-Ao. 
and so by (|9^) with i = at most 

5^S-(n-l)/2g(n-l)7rs 

for another constant 6„. So we have the estimate 

E^ 7 -Af-n/2 nTTS V^ x \ 

< aNOn s ' e 2^^ii ■ ■ ■ Ai„ 

where the smn on the right may be taken over all n-tuples (ii, . . . , i„). This sum is precisely equal to (triiT)" = s"^' 
Hence 






/2 nTTS 
2 



If we choose N > [n^ + n)/2 then we have 

2^2/ ^^^V2 



as desired. 
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C. Dyson's Continuum Model 

In ||l^ Dyson constructs a continuum Coulomb gas model |g^ for E/^{n; s). In this continuum model, 



where 



^ = ^ o / / Pix)p(y) log \^ - y\ d^dy 



2 
is the total energy, 

S — p[x) log p{x)dx 

is the entropy, p(x) — p{x) — 1 and p(x) is a continuum charge distribution on the line satisfying p{x) — > 1 as x — > ±00 
and p(x) > everywhere. The distribution p{x) is chosen to minimize the free energy subject to the condition 

p(x) dx = n . 

-s/2 

Analyzing his solution in the limit 1 << n « s, Dyson finds Ep{n\ s) ^ exp(— /3Wc) where 

Wc = -^- —(n + (5) + -n(n + J) + -n(n + 2<5) [log(— ) + -J (97) 

with (5=1/2- 1//3. 

We now compare these predictions of the continuum model with the exact results. First of all, this continuum 
prediction does not get the s~^l^ (for /3 = 2) or the s^^/* (for /3 = 1,4) present in all Ep{n;s) that come from 
the logTTs term in (|90|). Thus it is better to compare with the continuum prediction for rp{n; s). We find that the 
continuum model gives both the correct exponential behavior and the correct power of s for all three ensembles. 
Tracing Dyson's arguments shows that the power of s involving the n^ exponent is an energy effect and the power of 
s involving the n exponent is an entropy effect. Finally, the continuum model also makes a prediction (for large n) 
for the the Bp^n's (^2,n is given above). Here we find that the ratio of the exact result to the continuum model result 
is approximately n~^^^'^ ioi (3 — 2 and n^^^^^ for /S — 1,4. This prediction of the continuum model is better than it 
first appears when one considers that the constants themselves are of order n" /^ (/? = 2) and n" (/3 = 1, 4). 
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FIG. 1. Density of eigenvalues histogram for 25, 100 x 100 GOE matrices. Also plotted is the Wigner semicircle which is 
known to be the limiting distribution. 

FIG. 2. Density of eigenvalues histogram for 25, 100 x 100 symmetric matrices whose elements are uniformly distributed on 
[^1,1]. Also plotted is the Wigner semicircle distribution. 

FIG. 3. Level spacing histogram for 20, 100 x 100 GOE matrices. Also plotted is the Wigner surmise (g). 

FIG. 4. Level spacing histogram for 50, 100 x 100 symmetric matrices whose elements are uniformly distributed on [—1, 1]. 
Also plotted is the Wigner surmise (H). 

FIG. 5. Level spacing histogram for mixed data sets. Also plotted is the Poisson distribution. 
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